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Abstract 

Observations give strong support for the unification scheme of active galactic nuclei. 
The scheme is premised on toroidal obscuration of the central engine by dusty 
clouds that are individually very optically thick. These lectures summarize the torus 
properties, describe the handling and implications of its clumpy nature and present 
speculations about its dynamic origin. 
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1 Introduction 



The basic premise of the unification scheme is that every AGN is intrinsically 
the same object: an accreting supermassive black hole. This central engine 
is surrounded by a dusty toroidal structure so that the observed diversity 
simply reflects different viewing angles of an axisymmetric geometry. Since the 
torus provides anisotropic obscuration of the center, sources viewed face-on are 
recognized as "type 1" , those observed edge-on are "type 2" . The observational 
evidence for unification is covered in the lectures by C. Tadhunter and B. 
Peterson elsewhere in these proceedings. 

A scientific theory must make falsifiable predictions, and AGN unification does 
meet this criterion. Unification implies that for every class of type 1 objects 
there is a corresponding type 2 class, therefore the theory predicts that type 
2 QSO must exist. After many years of searching, the existence of QS02 has 
been established, thanks to the Sloan Digital Sky Survey, and their spectro- 
polarimetry even reveals the hidden type 1 nuclei at z as large as 0.6. This is a 
spectacular success of the unification approach. There are not that many cases 
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in astronomy — in fact, in all of science — where a new type of object has 
been predicted to exist and then actually discovered. In light of this success, 
it would be hard to question the basic validity of the unification approach. 
There is no reason, though, why the obscuring torus should be the same in 
every AGN; it is unrealistic to expect AGN's to differ only in their overall 
luminosity but be identical in all other aspects. 

The torus can be considered an acronym for Toro idal Obscuration Required 



by Unification Schemes. From basic considerations, iKrolik fc Begelmanl (119881 ) 
concluded that the obscuration is likely to consist of a large number of in- 
dividually very optically thick dusty clouds. However, the important roll of 
dumpiness was not always fully appreciated, for example, its effect on AGN 
classification statistics: type 1 and type 2 viewing is an angle-dependent prob- 
ability, not an absolute property. In these lectures I summarize the properties 
of the obscuring torus and try to speculate on how it might evolve with the 
AGN luminosity, i.e., its accretion rate. Since the main direct evidence for 
the torus comes from its IR emission, I start with a general discussion of dust 
emission (sec. 2). This is followed by a review of torus phenomenology (sec. 3), 
description of the handling of its dumpiness (sec. 4) and speculations about 
its dynamical origin (sec. 5). 



2 Dust— Generalities 



Dust grains are mixed with the gas in virtually all interstellar regions. The 
grains are small solid particles, built-up through collisions of many atoms (pri- 
marily Si, Mg, Al, C, O) sticking together. Grain interaction with radiation 
is comprised of two components: scattering and absorption] their sum is re- 
ferred to as extinction. Scattering occurs when radiation is simply reflected 
off the grain surface, while during absorption the incoming photon is actually 
absorbed by the dust grain. The interaction of dust with radiation follows the 
classical theory of scattering by small geometric structures characterized by a 
dielectric constant. The spectral variation of the cross-section a\ is primarily 
controlled by the relation between the radiation wavelength A and the grain 
size s. When A C s the grains block the radiation completely and a\ ap- 
proaches the constant diffraction limit. When A ^> s, ax decreases as a power 
law. Superimposed on this general behavior are spectral features reflecting in- 
ternal excitations of the grain material. Figure [1] shows the spectral behavior 
of each of the three cross sections for the "standard" composition of Galactic 
interstellar dust: 53% silicates and 47% graphite , with a grain size distribution 
follow ing a power-law from 0.005/im to 0.25/im (IMathis. Rumpl. fc Nordsiecki . 



19771 ) . Absorption is by far the dominant process for interaction of interstellar 



grains with radiation in the mid- and far-IR (wavelengths longer than ~ 2 
/^m); scattering can be safely ignored at these wavelengths. 
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Fig. 1. Spectral shape of interstellar dust extinction coefficient and its scattering 
and absorption components: owt = <7 sra , + o^bs- Standard i nters tellar grain mix 
with optical pr operties from Ossenkopf. Henning. &: Mathisl ( 19921 ) for the silicate 
component and Drainel ( 20031 ) for the graphite. 



The dust optical depth at wavelength A is r A = / n d a\d£, where <T\ is the appro- 
priate cross-section and is the dust number density. The mean projected 
area of dust grains shows the remarkably uniform correlation ridcry/un ~ 
5xlCT 22 cm 2 , where % is the hydrogen density and ay is the dust extinc- 
tion cross section at visual (0. 55/xm). This ratio is appro ximately the same 



in all Galactic directions (e.g.. ISparke fc Gallagher! . 120061 ). Dust obscuration 



with optical depth ry at visual therefore implies a gas column density Nn 
2xl0 21 r v cm" 2 . 



2.1 Dust Temperature 



The interaction of dust grains with radiation determines the state of their in- 
ternal excitation; collisional energy exchange with the gas is usually negligible 
in the dust energy balance (although it may be important for the gas). Dur- 
ing scattering the radiation is reflected off the grain surface and has no effect 
on the internal energy of the dust particle. However, when a grain absorbs 
an incoming photon it heats up and the energy is shared among the many 
internal excitation modes. In steady state, the grains must radiate away all 
the energy that they absorb. The balance between absorption and emission 
provides a condition of radiative equilibrium, which determines the grain tem- 
perature T: J a a hs\B\(T)d\ = J a^xJxdX, where B\ is the Planck function 
and J\ = J I\d\/4:7r is the angle-averaged intensity. The latter can be split 
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into its external and diffuse components according to J\ = J c \ + JdiffA- When 
the external radiation originates from a point source with luminosity L at a 
distance r, the bolometric magnitude of the external angle-averaged intensity 
on the illuminated surface of the dust distribution is J e = L/16irr 2 . Denote 
the Planck-average of the absorption coefficient cr a b S A by cr a p and its average 
with the spectral shape of the external radiation by a ac . The dust temperature 
equation is then 



o-a.p(T) 



L If 

°" bT4 = Tfi 2 + ^' Wllere $ = f / OabsA^difTAdA 



and where <7b is the Stefan-Boltzmann constant. The contribution of the dif- 
fuse radiation, contained in the term 5, couples this relation to the radiative 
transfer equation, and the dust temperature cannot be determined without 
solving the full dust radiative transfer problem. However, when the dust op- 
tical depth is small then the contribution of the diffuse radiation is negligible 
and S <C 1. In that case eq. [I] has an immediate solution, which gives roughly 
T ~ r ~ - 4 . Dust at temperature T emits predominantly around the peak of 
the Planck function, A ~ 10/xm (300 K/T). Because of the temperature de- 
cline with distance from the heat source, the wavelength of peak emission 
from distance r increases as ~ r 0A . As the dust optical depth increases, the 
temperature decline becomes steeper immediately behind the illuminated sur- 
face. However, once the temperature decreases to the point that the remain- 
ing dust is optically thin at wavelengths longer than the peak of the local 
Plan ck function, the optic ally-thin temperature variation is recovered (see, 
e.g., Ilvezic & Elitzurl . Il997h . 



As the grain size increases, the roughly flat portion of the cross section evident 
in figured] extends to longer wavelengths. The result is a larger cr a p and a lower 
grain temperature. Therefore, individual grains exposed to a heating radiation 
will have different temperatures, decreasing as the grain size increases, and the 
concept of a single "dust temperature" is meaningless. However, heating of all 
grains is coupled through the contribution of the diffuse radiation, and the 
different temperatures tend to equilibrate toward a single, common value in- 
side sources with large optical depths. A commonly used approximation is to 
replace the grain mixture with a single-type composite grain whose radiative 
constants are constructed from the mix average. This method reproduces ade- 
quately ful l calculations of grain mixtures , especially when o ptical depths are 



large (e.g.. lEfstathiou fc Rowan-Robinsonl . 1 19941 ; IWolfu2003l ). The handling of 



the dust optical properties is exact in this approach, the only approximation is 
in repla cing the tem peratures of the different grain components with a single 
average. IWolfi (120031 ) finds that the temperatures of different grains are within 
~ ± 10% of the temperature obtained in the mean grain approximation, and 
that these deviations disappear altogether when r v > 10. 
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Fig. 2. The dust temperature varies on the surface of an optically thick cloud in the 
radiation field of an external source from T max on the illuminated face to T m [ n on 
the dark side. Therefore the cloud's IR emission is anisotropic. The visible fraction 
of the illuminated surface area is determined by the position angle a with respect 
to the observer direction. 

Dust temperature distributions are profoundly different in clumpy and smooth 
environments. In smooth density distributions, dust temperature and distance 
from the radiative heating source are uniquely related to each other — given 
the distance, the dust temperature (either that of individual grain types or 
the common equilibrium value) is known, and vice versa. Clumpy media made 
up of discrete clouds that are individually optically thick and have dimensions 
that are much smaller than typical inter-cloud distances behave entirely dif- 
ferently. A cloud heated by a distant source will have a higher temperature on 
its illuminated face than on other surface areas (figure [2]) . The result is that 
in a clumpy medium, different dust temperatures coexist at the same distance 
from the radiation source, and the same dust temperature occurs at different 
distances — the dark side of a nearby cloud can be as warm as the bright side 
of a farther cloud. This has important implications for the IR emission from 
the AGN torus. 



Dust grains evaporate when they become too hot. Taking T su b = 1500 K as a 
typical sublimation temperature, the shortest wavelength of dust emission is 
~ 2/im; emission at shorter wavelengths requires dust hotter than ~1500 K, 
thus it necessarily reflects pure scattering. Inserting the spectral shape of the 
nuclear UV/optical radiation from active galaxies as the heat source in eq. [1] 
yields the AGN dust sublimation radius 

/ L \ 1/2 /1500K\ 26 
Rd — 0.4 { ) HF— ) P c ( 2 ) 



10 45 ergs -l J \ T t 



sub 



( INenkova et all l2008bl ) . This value of was determined from the tempera- 



ture on the illuminated face of an optically thick dusty cloud in the composite 
grain approximation. In fact, larger grains are cooler and can survive closer to 
the heat source. The sharp boundary denned in eq. [2] is an approximation. In 
reality, the transition between the dusty and dust-free environments is gradual 
because individual components of the mix sublimate at slightly different radii, 
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with the largest grains surviving closest to the AGN. 



2.2 The lOfim Feature 



Superimposed on the generally smooth spectral shape of the dust cross-section 
are two prominent features in the mid- infrared (fig. [1]). They are attributed 
to amorphous silicate grains that produce strong opacity peaks due to the 
Si-0 stretching and the O-Si-0 bending modes centered around 10 and 18 
fim, respectively. These silicate features, especially the 10/zm one, are visible in 
the observed spectral energy distribution (SED) of many astronomical objects, 
where they appear in either emission or absorption. This variety arises from 
radiative transfer effects as demonstrated in figure [31 which shows the results 
of exact model calculations^. 
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Fig. 3. Model results for dust with standard interstellar composition illuminated by 
black-body radiation with temperature of 10,000 K. The models labeled "sphere" 
have dust in a spherical shell with density profile 1/r 2 around a central point source. 
In those labeled "slab", parallel rays impinge on the dusty slab. The figure shows 
the normalized SED (J f\d\ = 1) of the radiation observed from outside the sphere 
and the non- illuminated face of the slab. The dust optical depth across the slab and 
along the shell radius is the same in each panel, as marked. 

1 The calculations were performed with the code DUSTY, which is publicly avail- 
able at |http: / /www.pa.uky.edu/~moshe/dusty| 
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The feature behavior is affected both by the dust optical thickness and its 
geometry. It can be understood with the aid of the simple solution of the ra- 
diative transfer equation for the emission from dust with uniform temperature 
T and optical depth r A : 



when the dust is exposed to external radiation with intensity J e , the transmis- 
sion term I e e~ Tx is added. This result shows that the emission from optically 
thin dust (r A <C 1) is proportional to the dust cross section, therefore the 
emergent radiation in this case follows the spectral shape shown in figure [TJ 
A superposition of different dust temperatures preserves the proportionality 
with T\, producing an emission feature in all optically thin sources. This fea- 
ture is evident in the top panel of figure [3J The optical depth of the solutions 
displayed in this panel is rather small (ry = 0.15), and the short wavelengths 
(A < 2/mi) emission is simply the input radiation, which is barely attenuated. 
Longer wavelengths show the dust emission, and the difference between the 
sphere and the slab results is an important reflection of the different tempera- 
ture profiles in the two geometries. In a spherical shell, the dust temperature 
decreases with radial distance because of the spatial dilution of the heating 
radiation (eq. [1]). The large range of temperatures spanned by the cool dust 
is reflected in the relatively broad shoulder of long-wavelength emission in 
the sphere solution. In contrast, the slab is exposed to parallel-rays radiation, 
which does not suffer any spatial dilution, therefore its temperature is uniform 
and the dust emission covers a more limited spectral range. 

Optically thin emission increases linearly with the overall optical depth. As 
this optical depth approaches unity, self-absorption sets in, higher order terms 
in the expansion of e~ TA become important and, as is evident from eq. [31 the 
emission saturates at the featureless Planck function: at constant tempera- 
ture, emission and self-absorption exactly balance each other, producing the 
black-body thermodynamic limit. Therefore, a single-temperature optically 
thick region can never produce a feature, neither in emission nor absorption. 
Furthermore, irrespective of optical depth, single-temperature dust will never 
produce an absorption feature; it can only produce an emission feature when 
optically thin. The emergence of the absorption feature visible in the bottom 
panel of figure [3] reflects the temperature stratification in actual dusty mate- 
rial: as the radiation propagates from hot regions toward the observer, it passes 
through cooler regions where it suffers absorption that is not balanced by the 
emission from these cooler regions. Thus, the temperature structure primarily 
determines the strength of the absorption feature at large optical depths. The 
models displayed in this panel have ry = 150. The input radiation is fully ex- 
tinguished within a short distance from the illuminated face, creating a large 





(3) 



when t\ 1; 
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temperature gradient close to the surface. Beyond that absorption layer, the 
temperature profile resembles the optically thin case and is controlled by the 
geometry — steep decline in the spherical shell and near constancy in the slab. 
This difference is reflected not only in the longer wavelength emission from 
the spherical shell but also the much greater depth of its 10/im absorption 
feature. 



This analysis shows that a temperature gradient is essential for an absorption 
feature; deep absorption requires dust geometry conducive to large gradients. 
Consider a cloud illuminated from outside by a radiative source. If the cloud's 
dimensions are much smaller than the distance to the source, then the heating 
flux is constant across the cloud's volume, just as in the slab model. Without 
radiative transfer effects, the dust temperature would be uniform throughout 
the cloud. In contrast, arranging the same dust in a geometrically thick shell 
around the same heating source produces a large temperature gradient because 
of the spatial dilution of the flux with radial distance. We can therefore expect 
the absorption feature to have only limited depth in the case of externally 
illuminated clouds. A deep feature requires the radiation source to be deeply 
embedded in dust that is thick both optically and geometrically. These results 
have important implications for the study of dust emission from active galaxies 



( ILevenson et al.l . 120071 ; ISirocky et all 120081 ) . 



3 Torus Phenomenology 



We start with a discussion of torus observations and the physical properties 
that can be inferred directly from them. 



3.1 IR Emission 



The primary evidence for the torus comes from spectropolarimetric observa- 
tions of type 2 sources, which reveal hidden type 1 emission via reflection off 
material situated above the torus opening. While compelling, this evidence is 
only indirect in that it involves obscuration, not direct emission by the torus 
itself. An obscuring dusty torus should reradiate in the IR the fraction of nu- 
clear luminosity it absorbs, providing direct evidence for its existence. Indeed, 
the continua from most AGNs show significant IR emission. The top panel of 
figure H] shows the composite type 1 spectra from a number of compilations. 
The optical/UV region shows the power law behavior expected from a hot disk 
emission. At A > 1/im, the SED shows the bump expected from dust emis- 
sion. The bottom panel shows the same data after subtracting a power law fit 
through the short wavelengths in a crude attempt to remove the direct AGN 
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Fig. 4. SED's of type 1 sources. Top: Average spectra from the indicated compi- 
lations; for the definition of Radio Quiet Quasars (RQQ) see Tadhunter's lectures. 
Bottom: The same SED's after subtracting a power law fit through the short wave- 
lengths (< 1/um). These AGN-subtracted SED's are similar to those observed in 
type 2 sources. 



component and mimic the SED from an equatorial viewing of these sources 
according to the unification scheme. Indeed, the AGN-subtracted SED's re- 
semble the observations of type 2 sources. Silicates reveal their presence in the 
dust through the 10 /xm feature. The feature appea rs generally in em ission in 
type 1 sources and in absorption in type 2 sources (IHao et all 120071 ). 



The torus dust emission has been resolved recently in 8-13 urn interferqmetr y 
with the VLTI. The first t orus was detected in N GC 1068 bv ljaffe et all fl2004h . 
the second in Circinus by lTristram et al.l (120071 ). The latter observations show 
evidence for the long anticipated clumpy structure. The dust temperature dis- 
tributions deduced from these observations indicated close proximity between 
hot (> 800 K) and much cooler (~ 200-300 K) dust. Such behavior is puzzling 
in the context of smooth-density calculations but is a natural consequence of 
clumpy models (see sec. 12. ip 
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Fig. 5. AGN classification according to unified schemes. Left: In a smooth-density 
torus, the viewing angle (from the axis) i = — a separates between type 1 and 
type 2 viewing. Right: In a clumpy, soft-edge torus, the probability for direct viewing 
of the AGN decreases away from the axis, but is always finite. 

3.2 Torus Obscuration 



There are clear indications that the optical depth in the torus equatorial plane 
is at least ry > 10. If the dust abundance in the torus is similar to Galactic in- 
terstellar regions, the equatorial column density is at least iV H > 2x 10 22 cm -2 . 

The classification of AGN into types 1 and 2 is based on the extent to which 
the nuclear region is visible, therefore source statistics can determine the an- 
gular extent of the torus obscuration. In its standard formulation, the uni- 
fication approach posits the viewing angle as the sole factor in determining 
the AGN type. This is indeed the case for a smooth-density torus that is 
optically thick within the angular width a (figure left sketch). All AGN 
viewed at < i < \n — o then appear as type 1 sources, while viewing at 
\ix — o < i < \k gives type 2 appearance. If f 2 denotes the fraction of type 
2 sources in the total population, then f 2 = sincr. This relation holds also 
when the column density declines smoothly with angle away from the equa- 
torial plane. The AGN is obscured from directions that have e _r 3> 1 and 
visible from those with e~ T <C 1. Because of the steep variation of e~ T with t, 
the transition between these two regions is sharp, occurring around the angle 
where r = 1. Denoting this angle a yields f 2 = sin a irrespective of the specific 
angular profile. Taking account of the torus dumpiness modifies this relation 
fundamentally, as is evident from the right sketch in figure [51 We discuss this 
modification below (sec. 14. 3ft . 



From statistics of Seyfert galaxies ISchmitt et al.l (120011 ) find that f 2 — 70%. 



Employing the standard f 2 = sincr relation th ey deduce a ~ 45°. The issue 



is currently unsettled because lHao et al.l (120051 ) have found recently that f 2 is 



only about 50%, or o ~ 30°. The fraction / 2 of obscured AGN decreases with 
bolometric luminosity. This has been verified in a large number of observations 
that estimate the luminosity dependence of either f 2 or fi (the fraction of 
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unobscured source s), or diffe r ences between t he luminosity functions of type 
1 and 2 AGN (see Irlao et all 120051 : ISimpsonl . 120051 : iMaiolino etaD . 120071 ). 



3.3 Torus Size 



Obscuration statistics provide an estimate of the torus angular width a. De- 
note by H the torus height at its outer radius R D , then H/R = tana; ob- 
scuration does not depend individually on either H or R D , only on their ra- 
tio. To determine an actual size one must rely on the torus emission. In the 
absence of high-resolution IR observations, early estimates of th e torus size 
came from theoretical analysis of the SED. iPier fc Krolikl (119921 ) preformed 
the first detailed calculations of dust radiative transfer in a toroidal geometry, 
and because of the difficulties in modeling a clumpy medium approximated 
the density distribution with a uniform one instead. They concluded that the 
torus has an outer radius R Q ~ 5-10 pc, but later speculated that this compact 
structure might b e embedded in a muc h larger, and more d i ffuse, torus with 
R Q ~ 30-100 pc (IPier & Krolikl . ll993T ). iGranato k Danese! (Il9941 ) extended 
the smooth-density calculations to more elaborate toroidal geometries. From 
comparisons of their model predictions with the observed IR emission at A 
~ 10-25 /im they concluded that the torus must have an outer radius R Q 
> 300-1000 p c , and that its radial density profile must be constant; later, 
Granato et all (119971 ) settled on hundreds of pc as their estimate for the torus 
size. Subsequently, R Q > 100 pc became common lore. 



The advent of high-resolution IR observations brought evidence in support of 
Pier fc Krolik's origina l proposal of compact torus dimensions (for details, see 
Nenkova et~aT1 . l2008bh . It can be argued that IR observations only determine 
the size of the corresponding emission region and that the actual torus size 
could in fact be much larger, but mid-IR flux considerations were the sole 
reason for introducing large sizes in the first place. It seems safe to conclude 
that at this time there is no compelling evidence that torus clouds beyond R Q 
~ 20-30.Rd need be considered. From eq. [2J a conservative upper bound on 
the torus outer radius is then R Q < 12 L\^ 2 pc, where L 45 = L/10 45 ergs" 1 . 



3.4 Torus Orientation and the Host Galaxy 



Toroidal obscuration is the cornerstone of the AGN unification scheme. Since 
the active nucleus is at the heart of a galaxy, might the obscuration not be 
attributed to the host galaxy? In fact, although galactic obscuration can af- 
fect individual sources, the strong orientation-dependent absorption cannot 
be generally attributed to the host galaxy because the AGN traced by 

the jet position angle, is randomly oriented with respect to the galactic disk 
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in Seyfer t galaxies (IKinnev et all |2000| ) an d the nuclear dust dis k in radio 



galaxies (ISchmitt et all 120021 ) . In addition, iGuainazzi et al.l (120011 ) find that 



heavily obscured AGN reside in a galactic environment that is as likely to be 
'dust-poor' as 'dust-rich'. 

NGC 1068, the archetype Seyfert 2 galaxy and one of the most studied active 
nuclei, is a revealing example. While the galaxy is oriented roughly face-on, 
the torus is edge-on to within ~ 5°. The nearly edge-on orientation of the torus 
and accreti on disk are indicated by the geometry and kinem atics of both wa- 
ter maser (iGreenhill Gwinnl. 119971; iGallimore et al.l . l200ll ) and narrow- line 
emission ( I Crenshaw fc Kraemerl . |2000| ). Some other interesting information 
about the AGN i mmediate surroundings in this source comes from molecular 
line observations. Schinnerer et al. fl2000h find from CO velocity dispersions 
that at R ~ 70 pc the height of the molecul a r clou d distribution is only H 
~ 9-10 pc, for H/R ~ 0.15. iGalliano et all (l2003h model the CO and H 2 



emission with a clumpy molecular disk with radius 140 pc and scale height 
20 pc, for the same H/R ~ 0.15. Thus, although resembling the putative 
torus, the distribution of these clouds does not meet the unification scheme 
requirement H/R ~ 1. Evidently, the detected m olecular clouds are lo cated 
in a thinner disk-like structure outside the torus. IGalliano et al.l (120031 ) find 
that this molecular disk is tilted roughly 15° from edge-on, much closer to the 
ori entation of the AGN t han of the host galaxy. Imaging polarimetry at 10/mi 
by iPackham et al.l (120071 ) shed some light on the continuity between the torus 
and the host galaxy's nuclear environments. The most recent high resolution 
adaptive optics observations of the nuclear region are presented by R. Davies 
elsewhere in these proceedings. 



3. 5 X-rays 



AGN are strong X-ray emitters, and the X-ray observations give overwhelm- 
ing evidence for the orientation-dependent absorption expected from the uni- 
fication scheme. The 2-10 keV X-ray continuum from Seyfert 1 nuclei is 
generally unattenuated. A small amount of neutral absorption affects some 
of these sources, and can be attributed to the host galaxy because it is 
usually characterize d by typical galactic column densities Nn< 10 21 cm -2 
( George et al. . 19981 ) J 2 ! In contrast, the continuum from Seyfe rt 2 galaxies is 



attenuated by obscuring columns with Nn ~ 10-10 cm (IBassani et al. 



19991 ; iRisaliti et allll999l ). In accordance with the unification scheme, the ab- 



sorption correct ed spectra and luminosities of type 2 sou rces are similar to 
those of type 1 ( Smith fc Done . 1996 ; Turner et al. , 1997bl ). From the correla- 



Spectral features give evidence for "warm absorbers" , clouds of ionized gas with 
2Vh ~ 10 21 -10 23 cm" 2 ; see George et al for details. 
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tions they find among various independent absorption indicators. iBassani et al. 
( 119991 ) conclude that the torus optical thickness along the line of sight is the 
parameter most closely correlated with the AGN X-ray properties. Promi- 
nent among these indicators is the Ka iron line, the strongest line in the 
4-10 keV X-ray spectrum o f AGN. In Seyfert 1 gala xies the line equivalent 
width is typically ~150 eV (INandra fc Pounds! . Il994l ). while in Seyfert 2 the 
equivalent widths are more br oadly distributed, ranging from about 100 eV 



to 1 keV (iTurner et all Il997al ) and even 5 keV in some Compton-thick AGN 



( ILevenson et all 120021 ) . The Seyfert 2 eq uivalent widths reflect processing in 



a tor oidal structure around the nucleus (jKrolik et all Il994t ILevenson et al. 



20021 ). 



Yet in spite of the overall correspondence between the optical and X-ray 
obscuration, there is a significant number of AGN for which the expected 
characteristics are different in the two bands. Although substantial X-ray ab- 
sorption is common among type 2 AGN, there are also unabsorbed X-ray 
sources that present only narrow emission lines in their optical spectra. Such 
cases can be explained wit h the observational se lection effect suggested by 



Severgnini et al.l (120031 ) and lSilverman et al.l (120051 ): in these sources, the opti- 



cal light of the host galaxy outshines the AGN continuu m and broad lines. This 
suggestion w a s supp orted by the subsequent studies of iPage et al.l (120061 ) and 



Garcet et al.l (120071 ) . In addition, some of the low- luminosity cases may reflect 



the disappearance of both torus obscuration and broad line emission (see sec. 
15.31 bellow). The opposite case, obscuration only in X-rays, exis ts too — there 
are type 1, broad l ine AGN with significant X-ray absorption (jPerola et al. 



2004 ; lEckart et al.l . 120061 ; I Garcet et all 120071 ) . Extreme cases include quasars 



whose optical spectrum shows little or no dust extinction w hile their X-ray con- 



tinuum is heavi l y affe cted by Compton thick absorption (IBraito et all 12004 



Gallagher etldl l2006j l. This cannot be attributed to observational selection 



effects. 



The differences observed in some sources between X-ray and UV/optical at- 
tenuation arise naturally from the different absorption properties of gas and 
dust. Dusty material absorbs continuum radiation both in the UV/optical and 
X-rays, therefore the dusty torus provides obscuration of both. But dust-free 
gas attenuates just the X-ray continuum, so clouds inside the dust sublima- 
tion radius will provide obscuration only in this band. Conclusive evidence 
for such obscuration comes from the short time scales for transit of X-ray ab- 
sorbing clouds across the line of sight, which es tablish the exi s tence of obscur- 
ing clouds inside the dust subli mation radius (jRisaliti et all 120021 ) . Extreme 
cases involve 4 hour variability (jElvis et al.l . 120041 ) and variations in absorbing 
column of more than 10 24 cm -2 within two days, indicating Compton thick 
(A^h > 10 24 cm~ 2 ) X-ray a bsorption from a single cloud in the broad-lines 



region (jRisaliti et all 120071 ) . These observations show that the torus extends 
inward beyond the dust sublimation point to some inner radius _R X < 
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Clouds at R x < r < partake in X-ray absorption but do not contribute 
appreciably to optical obscuration or IR emission because they are dust-free. 
Since every cloud that attenuates the optical continuum contributes also to 
X-ray obscuration but not the other way round, the X-ray absorbing column 
should be at least as lar ge as the UV/optical abs o rbing column, as observed 



(IMaccacaro et all Il982l ). Further, iMaiolino et al.l (120011 ) find that the X-ray 



absorbing column exceeds the reddening column in each member of an AGN 
sample by a factor ranging from ~3 up to ~100. Although this discrepancy 
could also arise from unusual dust-to-gas ratio or dust properties, the observa- 
tional evidence for X-ray absorption by clouds in the dust-free inner portion 
of the torus implies that in all likelih ood, these clouds proy ide the bulk of 
the absorption. This could explain the iGuainazzi et aJj (120051 ) finding that at 
least 50% of Seyfert 2 galaxies are Compton thick. 



3.6 X-ray vs IR 



IR flux measurements collect the emission from the entire torus area on the 
plane of the sky, originating from clouds along all rays through the torus. In 
contrast, X-ray attenuation is controlled by the clouds along just one partic- 
ular ray, the line of sight to the AGN. Therefore, even for the dusty portion 
of the obscuring column, the number of X-ray absorbing clouds can differ 
substantially from the torus average. Two type 2 sources with similar cloud 
properties and the same overall number of clouds following the same spatial 
distribution will have an identical IR appearance, yet their X-ray absorbing 
columns could still differ significantly. This can be expected to introduce a 
large scatter in torus X-ray properties among AGN with similar IR emission. 
It may help explain why the SEDs show only moderate variations in the in- 
frared that are n ot well correlated with the X-ray absorbing columns (e.g., 



Silva et all 120041 ) . 



4 Clumpiness — Handling and Implications 



The clumpy nature of the environment around the central black-hole is crucial 
for understanding various aspects of A GN observations. We start by present- 
ing a general formalism developed in iNenkova et al.l (120021 . l2008al ) for han- 
dling clumpy media. The original development, described here, applies only 
to continuum emission and takes all clouds to be identical, for simplicity. The 
for malism was extended to line emi s sion a nd to a mixture of cloud properties 
by IConwav. Elitzur. fc Parral (120051 . l2008h . 
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4-1 Clumpy Radiative Transfer 



Fig. 6. A region populated by clouds of size R c . Along ray s, the photon mean-free- 
path is £ and the number of clouds per unit length is Nq = £ ■ 

Figure [6] sketches a region where the matter is concentrated in clouds. In- 
dividual clouds are characterized by their size R c , the cloud distribution is 
specified by the number of clouds per unit volume no It is useful to introduce 
the number of clouds per unit length Nq = ricA c = l~ x where A c is the cloud 
cross-sectional area and I is the photon mean free path for travel between 
clouds. Denote by V c the volume of a single cloud and by <fi the volume filling 
factor of all clouds, i.e., the fraction of the overall volume that they occupy. 
The medium is clumpy whenever 

4> = n c K « 1 ■ (4) 



In contrast, the matter distribution is smooth, or continuous, when (f) ~ 1. 
Since V c ~ A C R C , the dumpiness condition is equivalent to 

(f) = N C R C < 1, or R c < 1. (5) 



The dumpiness criterion is met when the mean free path between clouds 
greatly exceeds the cloud size. Under these circumstances, each cloud can 
be considered a "mega-particle", a point source of intensity S C) \ and optical 
depth T\. The intensity at an arbitrary point s along a given path can then be 
calculated by applying the formal solution of radiative transfer to the clumpy 
medium. The intensity generated in segment ds' around a previous point s' 
along the path is S c ,\(s')Nc(s')ds' . Denote by J\f(s',s) = J*, Neds the mean 
number of clouds between s' and s and by P esc (s',s) the probability that 
the radiation from s' will reach s without absorption by intervening clouds. 
The number distribution of clouds between s' and s follows Poisson statistics 



around the mean J\f(s', s), and lNatta &: Panagial (119841 ) show that 



P esc (s',s) = e-^ s '> s \ where t x (s', s) = M(s', s)(l - e^). (6) 
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The intuitive meaning of this result is straightforward in two limiting cases. 
When T\ < 1 we have t\(s', s) ~ Af(s f , s)t x , which is the overall optical depth 
between points s' and s; that is, dumpiness is irrelevant when individual 
clouds are optically thin, and the region can be handled with the smooth- 
density approach. It is important to note that Mt\ can be large — the only 
requirement for this limit is that each cloud be optically thin. The opposite 
limit T\ > 1 gives P esc (s', s) ~ e --^0'> s ). Even though each cloud is optically 
thick, a photon can still travel between two points along the path if it avoids 
all the clouds in between. With this result, the intensity at s generated by 
clouds along the given ray is 



This relation is the exact analog of the formal solution of standard radiative 
transfer in continuous media, to which it reverts when the cloud sizes are 
decreased to the point that they become microscopic particles; in that case 
t <C 1 for each particle and i~ x is the standard absorption coefficient. The 
result does not involve 0; only Nq enters. A complete formalism that would 
not invoke the assumption 1 would lead to a series expansion in powers 
of 0, and our derived expression would yield the zeroth order term in that 
expansion. In fact, detailed Monte Carlo simulations show that, to within a 
few percents, our result describes adequately clumpy media with as large as 
10%, as shown in figure [3 Since the intensity calculations are independent of 
the volume filling factor, modeling results do not provide any information on 
this quantity, nor do they provide separate information on either R c or nc, 
only on Nq- In complete analogy, the radiative transfer problem for smooth 
density distributions does not involve separately the size of the dust grains 
or their volume density, only the combination ndCd, which determines the 
absorption coefficient and which is equivalent to Nq- 

An intensity calculation requires the cloud source function S c> \ at every loca- 
tion in the clump distribution (eq. [7j). Since the clump emission is anisotropic, 
as is evident from figure [2J S C}X depends not only on location but also on the 
observer's position angle a. It is convenient to split the population of clumps 
to those that are directly exposed to the AGN radiation, with a source func- 
tion S^ A , and those shadowed by intervening clouds. The latter are heated 
only indirectly by the emission from all other clouds and their source function 
is denoted S^. A . At location (r, where j3 is angle from the equatorial plane, 
the mean number of clouds to the AGN is N(r,0) = JqNc (r, [3) dr and the 
probability for unhindered view of the AGN is p(r, p) = e~ M ^\ The general 
expression for the cloud overall source function is thus 




(7) 



S c , x (r, a, (3) = p(r, f3)S* x (r, a) + [1 - p(r, (3)} S^r, a, (3) 
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Fig. 7. Averaged transmission, exp(— r av ), through a clumpy region with filling factor 
(f>. The optical depth of a single cloud is to, as marked on top each panel; all To > 10 
are equivalent. Horizontal lines show the theoretical small-</> limit (eq. [6]); they are 
labeled by the mean number of clouds along a line through the volume. The results 
of Monte Carlo simulations are shown with filled circles whose sizes are comparable 
to the numerics error bars. Gray curves trace the loci of constant cloud sizes. The 
analytical expression adequately describes clumpy regions with 4> as large as 10% 
(jConwav et all l2008h . 



The source functions S^ x can be calculated from a straightforward solution of 
the standard radiative transfer problem once a geometrical shape is assumed 
for the individual clouds, a calculation that is independent of the cloud distri- 
bution. The indirect heating of shadowed clouds, on the other hand, involves 
the radiation from all other clouds and thus cannot be calculated without 
solving the full problem. This is the difficulty presented by diffuse radiation 
in standard radiative transfer. Since the diffuse heating is dominated by the 
clouds that are directly heated by the AGN, the S^ x source functions can 
serve as the starting point for an iterative solution of the full problem: From 
the emission of the directly heated clouds devise a first approximation for the 
diffuse radiation field. Next, place clouds in this radiation field and calculate 
their emission to derive a first approximation for the source functions of in- 
directly illuminated clouds and, from eq. the composite source function at 
every location. In successive iterations, add to the AGN direct field the cloud 
radiation calculated from e q. [7J and repeat the process until convergence. 



Nenkova et al.l (120021 . l2008al ) performed only the first step of this procedure, 



an approximation that leaves out the feedback effect of the clouds on their 
own radiative heating. Carrying out higher orders of the iteration scheme is 
an open challenge. 
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Fig. 8. Model spectra for dust IR radiation of a clumpy torus. In the bottom panel 
the cloud distribution follows eq. [9] with Mo = 5, a = 45° and q = 1 in the radial 
range < r < 30Rd- Each cloud has an optical depth ry = 60. Models in the top 
panel are the same, except that the angular distribution is uniform within \/3\ < a. 
Different curves show viewing angles that vary in 10° steps f rom pole-on (i = 0°) to 
edge-on (i = 90°). Fluxes are scaled with F A gn = L/4ttD 2 . (jNenkova et all l2008bl ) 



4-2 Clumpy Torus Emission 



The only distribution required for calculating the intensity from eq. [7|is iVc, 
the number of clouds per unit length. The torus axial symmetry implies that 
Nq depends only on the distance r from the center and the angle (3 from the 
equatorial plane (the complementary of the standard polar angle). Denote by 
Ao the mean of the total number of clouds along radial rays in the equatorial 
plane. We parametrize the angular profile of the cloud distribution as a Gaus- 
sian with width parameter a and the radial distribution as a declining power 
law with index q, so that 

iVc(r,/?) = £Ar e-^ 2 (^J (9) 
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Fig. 9. Fits of clumpy torus models to small-aperture (0.4") Gemini 8-13 //m spec- 
trum of NGC 1068, shown with heavy solid line. Symbols and thin solid line between 
8-13 //m mark all other data, which were not utilized in the fits. Filled symbols 
denote measurements with aperture a < 0.5", open symbols a > 0.5" . Model spec- 



tra a re for equatorial viewing and various parameters (for details see iMason et al. 
20061 ) . 



where C = (/ 1 y dy/y 11 )' 1 is a dimensionless constant, ensuring the normaliza- 
tion Ao = / Nc(r, (3 = 0)dr. The Gaussian can be taken as representative of 
soft-edge angular distributions, illustrated on the right in figure [5j Sharp-edge 
angular distributions can be parametrized with a step function with cutoff 
at (3 — a. Figure M shows the results of model calculations for these two ge- 
ometries with the same set of parameters. The sharp-edge geometry produces 
a bimodal distribution of spectral shapes, with little dependence on view- 
ing angle other than the abrupt change that occurs between the torus open- 
i ng and the obscured region . This SED dichotomy conflicts with observations 



(jAlonso-Herrero et all 120031 ). In contrast, the soft-edge Gaussian distribution 



produces a larger variety in model spectral shapes, wit h a smooth, continu- 



ous de pendence on i, in agreement with the findings of lAlonso-Herrero et al. 



(120031 ). It is gratifying that the more plausible soft-edge distribution also pro- 



duces better agreement with observations. 

Clumpy torus modeling was employed in the analysis of spatially-resolved, 
ne ar-diffraction- l imite d 10 /an spectra of the nucleus of NGC 1068, obtained 



by lMason et al.l (120061 ) with Michelle, the mid-IR imager and spectrometer on 
the Gemini North telescope. Mason et al first compared the models with the 
8-13 /im Michelle spectrum of the central 0.4" core to determine the best- 
fit parameters. The model predictions outside this spectral range were then 
compared with all other observations, and figure M shows the results of this 
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comparison. As noted above, the torus in this AGN is nearly edge-on (sec. 
I3.4p . and this is the first model to reproduce the observed near-IR flux with 
such orientation; in contrast, smooth-density models were forced to assume, 
contrary to observations, viewing from 22°-25° above the equ atorial plane in 



order to bring into view the warm face of the torus backside (IGranato et al. 



19971 ; iGratadour et all 120031 ; iFritz et all 120061 ) . It is noteworthy that all flux 



measurements with apertures < 0.5" are in good agreement with the model 
results. However, the flux collected with larger apertures greatly exceeds the 
model predictions at wavelengths longer than ~4/zm. This discrepancy can be 
attributed to IR emission from nearby dust outside the torus — Mason et al 
show that the torus contributes less than 30% of the 10 /im flux collected with 
apertures > 1". The bulk of the large-aperture flux comes at these wavelengths 
from dust in the ionization conegjj; while less bright than the torus dust, it 
occupies a much larger volume. On the other hand, the torus dominates the 
emission at short wavelengths; at 2 /im, more than 80% of the flux measured 
with apertu res > 1" comes from the torus even though its image size is less 



than 0.04" (j Wei gelt et all 12004 ) . This highlights a difficult problem that af- 
flicts all IR studies of AGN. The torus emission can be expected to dominate 
the observed flux at near IR because such emission requires hot dust that ex- 
ists only close to the center. But longer wavelengths originate from cooler dust, 
and the torus contribution can be overwhelmed by the surrounding regions. 
There are no easy solutions to this problem. 

All in all, clumpy torus models seem to produce SED's that are in reason- 
able overall agreement with o bservations for the following range of parameters 
flNenkova et all l200l l2008al lbh : 



• Number of clouds, on average, along radial equatorial rays Ao = 5-15 

• Gaussian angular distribution with width a = 30°-50° 

• Index of power law radial distribution q = 1-2 

• Optical depth, at visual, of each cloud Ty = 30-100 

1 11 

• The torus extends from dust sublimation at Ra = 0.4L 4 5 pc to an outer 
radius R Q > 5Rd 

The silicate features display in clumpy models a complex behavior that can- 
not be reproduced with smooth-density distributions. Figure [10] shows some 
model results for the spectral region containing the features. The depth of 
the absorption features is never large in spite of the huge optical depths of 
some of the displayed models — a torus with Ty = 150 and Ao = 10 has an 
equatorial optical depth of 1500 at visual. The relatively shallow absorption 
reflects the small temperature gradient across an externally illuminated cloud, 
similar to the slab case (see sec. 12.21 fig. [3]). Indeed, a striking characteristic 
of all AGN spectra is the absence of any deep 10 /im absorption features in 



For discussion of ionization cones see lectures by Netzer, Peterson and Tadhunter. 
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Fig. 10. Spectral shape of the silicate 10 and 18 /jm features: Fx is the torus emission 
in the 5-30 fim region and F c- \ is the smooth underlying continuum obtained by a 
spline connecting the feature-free segments in this spectral region. All models have 
q = 2, R Q = 30i?d 5 o = 45°, and Mo as marked in each panel. Curves correspond 
to different ry, as labeled. Panels on the left correspond to pole-on viewing, on 
the right to edge-on; n ote the different scales of the vertical axes in the two cases. 
(|Nenkova et all l2008bh 



contrast with Ul traluminous IR G alaxies (ULIRGs) where the features reach 
extreme depths (IHao et all 120071 ). This different behavior in ULIRGs can be 
attributed to deep embedding in a dust dis tribution that is smooth, rather 



than clumpy (ILevenson et al.l 120071 ; see also ISpoon et al.l 120071 . ISirocky et al 
20081 ) . A most peculiar result is the emergence of the 10/zm feature in emission 
at equatorial viewing when ry is large for all Ao, and when J\f is small for all 
Ty. This could e xplain intrigu i ng Sp itzer observations of seven high-luminosity 
type 2 QSOs by lSturm et al.l (120061 ). Although the individual spectra appear 
featureless, the sample average spectrum shows the 10 mic feature in emission. 
Therefore, if this finding is verified it could indicate that the optical depths of 
torus clouds perhaps are larger in QSOs than in Seyfert galaxies or that the 
cloud number Ao decreases as the luminosity increases. 



While the parameters listed above lead to SED's compatible with observations, 
some additional considerations can further restrict the acceptable range. For 
example, in spite of the huge difference in AGN obscuration between type 
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1 and 2 sources, comparisons of th eir IR emission seem to indicate a sur- 



prisingly low amo unt of anisotropy (ILutz et all 12004 ; iBuchanan et all 12006 



Horst et all 120061 ). Models with q = 2 tend to produce less variation with 
viewing- angle than those with q = 1. Therefore, observations that will place 
reliable tight constraints on the IR emission anisotropy may indicate that q = 
2 provides a more appropriate description of the radial distribution than q = 
1. Finally, all model calculations were performed with standard Galactic ISM 
dust, which seems to provide satisfactory results. Current data do not provide 
any compelling reason for drastic changes in the dust composition. 



4-3 Clumpy Unification 



Because of dumpiness, the difference between types 1 and 2 is not truly an is- 
sue of orientation but of probability for direct view of the AGN (figure right 
sketch); AGN type is a viewing- dependent probability. Since that probability is 
always finite, type 1 sources can be detected from what are typically consid- 
ered type 2 orientations, even through the torus equatorial plane: if Mq = 5, for 
example, the probability for that is e -5 = 1/148 on ave rage. This might offer 



an ex planation for the few Seyfert galaxies reported by lAlonso-Herrero et al. 



( 120031 ) to show type 1 optical line spectra together with 0.4-16 jum SED that 
resemble type 2. Conversely, if a cloud happens to obscure the AGN from an 
observer, that object would be classified as type 2 irrespective of the viewing 
angle. In cases of such single cloud obscuration, on occasion the cloud may 
move out of the line-of-sight, creating a clear path to the nucleus and a tran- 
sition to type 1 spectrum. Such transitions betw een type 1 and t ype 2 line 



spectra have been observed in a few sources (see lAretxaga et all Il999l . and 
references therein). It is worth while to conduct monitoring observations in an 
attempt to detect additional such transitions. The most promising candidates 
would be obscured systems with relatively small X-ray obscuring columns, 
small torus size s (low A GN luminosity; see eq. [2]) and large black-hole masses 



(INenkova et alll2008ri ). 



A sharp-edge clumpy torus has f% = (1 — e - ^ ) sin a, which is practically in- 
distinguishable from the smooth-density relation when Mq exceeds ~ 3-4 (see 
sec. 13. 2j) . However, the situation changes fundamentally for soft-edge clumpy 
distributions because at every viewing angle, the probability of obscuration 
increases with the number of clouds. As is evident from figure [TT| the Gaussian 
distribution produces a strong dependence on Mq and significant differences 
from the sharp-edge case. Since the sharp-edge angular distribution is ruled 
out by observations, the fraction of obscured sources depends not only on the 
torus angular width but also on the average number of clouds along radial rays. 
While the fraction $2 = 70% requires a = 45° in the smooth-density case, it 
implies a = 33° in a Gaussian clumpy torus with Mq = 5 clouds. 
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Fig. 11. AGN statistics: The fraction f 2 of obscured sources as a function of the torus 
width parameter a. In torii with either smooth density distribution or a sharp-edge 
clumpy one, this fraction is determined uniquely by a, and is shown with the dashed 
line. In contrast, in a clumpy torus with Gaussian angular distribution (eq. [9]), f% 
depe nds on both a and the cloud number Mo, which is marked on the various solid 
lines dNenkova et all l2008bh . 



Observations show that f 2 decreases with L. This has been interpreted as sup- 
port for the "receding torus" model, in which a decreases with L (Simpson 
2005 and references therein). However, all the quantitative analyses performed 
thus far for the Independence of f 2 were based on the sharp-edge expression 
f 2 = sin a. Removing this assumption affects profoundly the foundation of the 
receding torus model because the dependence on the number of clouds neces- 
sitates analysis with two free parameters, therefore a cannot be determined 
without Mo- A decrease of Mo with L at constant o will also produce a decrease 
in f 2 , the same effect as a decrease of a (figure ITT]) . An observed trend of f 2 
with L may arise from a dependence on either a or Mo or both. There is no 
obvious a-priori means for deciding between the various possibilities. 



4-4 Torus Mass; Total Number of Clouds 



The total mass in torus clouds can be written as M torus = m p Nji JNc(r, (3) dV, 
where m p is the proton mass and N-r is the column density of a single cloud; 
note that M torus does not involve the volume filling factor. Denote Y = R Q / 
the torus radial thickness. With the cloud distribution from eq. [9] and taking 
for simplicity a sharp-edge angular distribution, so that the integration is 
analytic, the torus mass is M tor us = 47rm p sin a N tOTUS Rd 2 Y I q (Y), where I q = 
1, K/(21nY) and hY for q = 2, 1 and 0, respectively. Taking R d from eq. [21 
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the mass ratio of the torus and the central black hole is 

M^p = 2 xl0- 4 -^- sin aN torns<23 YI q (10) 



where Z/Edd is the Eddington luminosity and A^ orus 2 3 is the equatorial column 
density in 10 23 cm -2 . Since the Eddington ratio L/L^a is always < 1 and 
the radial thickness Y is in all likelihood no larger than ~20-30 (sec. 13 . 3[) . 
the torus mass is always negligible in comparison with M. when q — 2. If the 
radial cloud distribution is flatter, eq. [10] may constrain the torus properties 
to keep its mass below that of the black-hole. 

The total number of torus clouds is n tot = J n c dV = J(N C / A c )dV. This is the 
only quantity of interest that depends explicitly on the cloud size. Equivalently, 
R c can be replaced by the volume filling factor 0, since inserting eq. [9] into eq. [5] 
yields R c = C^R^/Mq at the torus inner edge. If is constant throughout the 
torus then n tot ~ Af\/(f) 2 for the 1/r 2 distribution, independent of the torus 
radial thickness Y. For example, if the volume filling factor is 10%, in order to 
encounter J\f = 5-10 clouds along each radial equatorial ray, the torus must 
contain n to t — 10 4 -10 5 clouds, 



5 What is the Torus? 



In the ubiquitous sketch by lUrry fc Padovanil (119951 ). shown numerous times 
also in this summer school, the AGN central region, comprised of the super- 
massive black hole (SBH), its accretion disk and the broad-line emitting 
clouds, is surrounded by a large doughnut-like structure — the torus. This hy- 
drostatic object is a separate entity, presumably populated by molecular clouds 
accreted from the galaxy. Gravity controls the orbital motions of the clouds, 
but the origin of vertical motions capable of sustaining the "doughnut" as a 
hydrostatic structure whose height is comparable to its radius was rec o gnize d 
as a problem since the first theoretical study by iKrolik fc Begelmanl (119881 ). 
An entir ely different scenario, hy drodynamic rather than hydrostatic, was pro- 
posed by lEmmering et al.l (119921 ). It involves the outflow of clouds embedded in 
a hydromagnetic disk wind, avoiding the vertical support problem. The clouds 
are accelerated by the system rotation alo ng magnetic field lines anc hored in 
the disk in a manner first described in the iBlandford fc Paynd (119821 ) seminal 
work. In this approach the torus is merely a region in the wind which happens 
to provide the required toroidal obscuration because the clouds there are dusty 
and optically thick. That is, instead of a hydrostatic "doughnut", the torus 
is just one region in a clumpy wind coming off the black-hole accretion disk. 
The wind paradigm wa s strongly advocated by A. Konigl and his associates 



(IKonigl fc Kartid . I1994J : IKartie fe Konigll . Il996l : IKartie et all Il999l : Everett 
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20051 ) and by iBottorff et al.l (119971 . |2000| ) but did not gain wide acceptance 
until recently because of the general perception that the torus must extend to 
more than ~100 pc, a distance too large to support the required disk wind. 



Recent observations suggest that the disk wind scena r io mi ght offer the bet- 
ter paradigm for the AGN torus (jElitzur Shlosmanl . 120061 ). High- resolution 
IR observations indicate that the torus size might actually be no more than 
a few pc (sec. 13.31) . This compact size places the torus clouds well inside 
the black hole sphere of influence, where the SBH gravity dominates over 
the galactic bulge. The outer boundary of this region is where the SBH and 
the bulge induce gravitational rotations with equal angular velocities, i.e., 
the radius within which the bulge mass is equal to that of the black hole. 
Consider a spherical bulge that induces a linearly rising rotation curve wit h 



Q ~ lkms 1 pc as is typical of AGN host galaxies (ISofue et al.l . 119991 ) 



and an SBH with a mass of M.7XlO 7 M at its center. The SBH will dom- 
inate the gravitational motions within a radius 35pcx(M.7 / /f2 : j ! ) 1 / 3 , where 
Qi —Q/(l kms" 1 pc -1 ). Since the torus is well within this sphere of influence, 
its dynamic origin is determined in all likelihood not by the host galaxy but by 
the central engine and its accretion disk. There is mounting evidence for cloud 
outflow in these regions with the geome try and kinematics of disk winds (e.g., 
Elvisl . I2004J ; iGallagher fc Everettl . 120071 ) , in accord with the outflow paradigm. 



Two different types of observations show that the torus is a smooth continu- 
ation of the bro ad lines region (BLR), n ot a separate entity. IR reverberation 
observations by iSuganuma et al.l (120061 ) measure the time lag of the dust ra- 
diative response to temporal variations of the AGN luminosity, determining 
the torus innermost radius. Their results show that this radius scales with lu- 
minosity as L 1 / 2 and is uncorrelated with the black hole mass, demonstrating 
that the torus inner boundary is controlled by dust sublimation (eq. [2]), not by 
dynamical processes. Moreover, in each AGN for which both data exist, the 
IR time lag is the upper bound on all time lags measured in the broad lines, a 
relation verified over a range of 10 6 in luminosity. This finding shows that the 
BLR extends outward all the way t o the inner boundary of the dusty torus, 
validating the iNetzer &: Laorl (119931 ) proposal that the BLR size is bounded 
by dust sublimation. The other evidence is the finding by lRisaliti et al.l (120021 ) 
that the X-ray absorbing columns in Seyfert 2 galaxies display time varia- 
tions caused by cloud transit across the line of sight. Most variations come 
from clouds that are dust free because of their proximity (< 0.1 pc) to the 
AGN, but some involve dusty clouds at a few pc. Other than the different time 
scales for variability, there is no discernible difference between the dust-free 
and dusty X-ray absorbing clouds, nor are there any gaps in the distribution. 
These observations show that the X-ray absorption, broad line emission and 
dust obscuration and reprocessing are produced by a single, continuous distri- 
bution of clouds. The different radiative signatures merely reflect the change in 
cloud composition across the dust sublimation radius R^. The inner clouds are 
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dust free. Their gas is directly exposed to the AGN ionizing continuum, there- 
fore it is atomic and ionized, producing the broad emission lines. The outer 
clouds are dusty, therefore their gas is shielded from the ionizing radiation, and 
the atomic line emission is quenched. Instead, these clouds are molecular and 
dusty, obscuring the optical/UV emission from the inner regions and emitting 
IR. Thus the BLR occupies r < while the torus is simply the r > R& region. 
Both regions absorb X-rays, but because most of the clouds along each radial 
ray reside in its BLR segment, that is where the bulk of the X-ray obscuration 
is produced. Since the unification torus is just the outer portion of the cloud 
distribution and not an independent structure, it is appropriate to rename it 
the TOR for Toroidal Obscuration Region. The close proximity of BLR and 
TOR clouds should result in cases of partial obscuration, possibly leading to 
observational constraints on cloud sizes. 



The merger of the ionized and the dusty clouds into a single population fits 
naturally into the TOR outflow paradigm. The AGN accretion disk appears to 
be fed b y a midplane influx of c old, clumpy material from the main body of the 



galaxy (IShlosman et al.l . Il990l . and references therein) . Approaching the cen- 
ter, conditions for developing hydromagnetically- or radiatively-driven winds 
above this equatorial inflow become more favorable. The disk-wind rotating 
geom etry provides a natu r al cha nnel for angular momentum outflow from the 
disk (IBlandford fc Paynd . Il982l ) and is found on many spatial scales, from 
protostars to AGN. The composition along each streamline reflects the origin 
of the outflow material at the disk surface. The disk outer regions are dusty 
and m olecular, as observed in water masers in some edge-on cases ( iGreenhilll . 



20051 ). and clouds uplifted from these outer regions feed the T OR. Such clouds 



have been detected in water maser observation s of Circinus (IGreenhill et al. 



20031) and NGC 3079 flKondratko et al.l . |2005|). The Circinus Seyfert 2 core 
provides the best glimpse of the AGN dusty /molec ular component. Wate r 
masers trace both a Keplerian disk and a disk outflow (IGreenhill et al.l . 120031 ) . 
Dust emission at 8-13/im shows a disk embedded in a slightly cooler and 
larger, geometrically thick torus (ITristram et al.l . 120071 ) . The dusty disk coin- 
cides with the maser disk in both orientation and size. The outflow masers 
trace only parts of the torus. The lack of full coverage can be attributed to 
the selectivity of maser operation — strong emission requires both pump action 
to invert the maser molecules in individual clouds and coincide nce along the 



line o f sight in both position and velocity of two maser clouds ( IKartje et al. 



19991 ). In NGC 3079, four maser features were found significantly out of the 



plane of the maser-traced disk yet their line -of-sight velocities reflect the veloc- 
ity of the most proximate side of the disk. iKondratko et al.l ( 20051 ) note that 



this can be explained if, as proposed by IKartje et al.l (I1999I ). maser clouds 



rise to high latitudes above the rotating structure while still maintaining, to 
some degree, the rotational velocity imprinted by the parent disk. Because 
the detected maser emission involves cloud-cloud amplification that requires 
precise alignment in both position and velocity along the line-of-sight, the 
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discovery of four high-latitude maser features implies the existence of many 
more such clouds partaking in the outflow in this source. Moving inward from 
the dusty/molecular regions, at some smaller radius the dust is destroyed and 
the disk composition switches to atomic and ionized , prod ucing a double-peak 



signature in some emission line profiles (jEracleousl . 12004 ) . The outflow from 



the atomic/ionized inner region feeds the BLR and produce s many atomic lin e 



signatures, including evidence for the disk wind geometry (lHall et al.1 . 120031 ) 



The outflow paradigm for the obscuring torus provides what can be called a 
"Grand Unification Theory" for AGN. This picture requires only the accretion 
disk and the clumpy outflow it generates — there is no need to invoke any 
additional, separate structure such as a hydrostatic "doughnut". Cloud radial 
distance from the AGN center and vertical height above the accretion disk 
explain the rich variety of observed radiative phenomena. In both the inner 
and outer outflow regions, as the clouds rise and move away from the disk 
they expand and lose their column density, limiting the vertical scope of X-ray 
absorption, broad line emission and dust obscuration and emission. The result 
is a toroidal geometry for both the BLR and the TOR. With further rise and 
expansion, the density decreases so the ionization parameter increases, turning 
the clouds into members of the warm absorber population (see lectures by H. 
Netzer elsewhere in these proceedings). Similar considerations can be invoked 
to explain the broad absor ption lines observed in some quasars (BAL/QSO; 



Gallagher fc Everettl . 120071 ). We proceed now to discuss some consequences of 



the TOR outflow scenario. 



5. 1 TOR Cloud Properties 



The only property of individual clouds constrained from the IR modeling is 
their optical depth. It should lie in the range Ty ~ 30-100 (sec. 14.21) . i.e., the 
cloud column density is N-r ~ 10 22 -10 23 cm -2 assuming standard dust-to-gas 
ratio. Clouds uplifted into the wind expand while moving away from the disk 
and their column density decreases. A cloud starting with A^h ~ 10 23 cm -2 
ceases to partake in obscuration when its column is reduced by ~ 100 so that 
its Ty drops below unity. Clouds starting with a smaller column will rise to a 
smaller height before losing their obscuration and leaving the TOR popula- 
tion. The obscuration properties constrain only the cloud column density, the 
product of the cloud's size and density. The latter can be constrained sepa- 
rately when we note that the cloud must be able to withstand the black-hole 
tidal shearing effect. Consider the tidal torque of the central SBH at a distance 
r pc in pc, where the Keplerian period is £k = 3x 10 4 M. 7 1/ ' 2 rp/ 2 yr. To prevent 
a cloud with density p = m p n from shearing, it must be at least partially 
confined by its own gravity and/or the ambient magnetic field B. In the for- 
mer case, the characteristic Jeans timescale is tj ~ (Gp)^ 1 ^ 2 = SxlO 4 ?^ 1 ^ 2 yr, 
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where n 7 = nj (10 7 cm -3 ). Therefore, resistance to tidal shearing (tj < £k) sets 
a lower limit on the density of the cloud, leading to upper limits on its size R c 
and mass M c \\ 

nr>^f, R c < 10 16 ^icm, M d < 7x10-^,23^6 M (11) 



where R cW = _R c /(10 16 cm) and A^ H ,23 = A^" H /(10 23 cm -2 ). The resistance to 
tidal shearing does not guarantee confinement against the dispersive force 
of internal pressure. Self-gravity cannot confine these clouds because their 
low masses provide gravitational confinement only against internal motions 
with velocities < 0.1 kins" 1 . A corollary is that these clouds cannot collapse 
gravitationally to form stars. While self-gravity cannot hold a cloud together 

1/2 

against dispersal, an external magnetic field B ~ 1.5 05% mG would suffice 
if the internal velocity dispersion is c^xlkms -1 . Clouds with these very same 
properties can explain the masers detected in Circinus and NGC 3079, adding 
support to the suggestion that the outflow water masers are yet another man- 
ifestation of the dusty, molecular clouds that make up the torus region of the 
disk-wind. Densities of H2O masers are frequently quoted as ~ 10 8 -10 9 cm -3 , 
but these are the optimal densities to produce the highest possible bright- 
ness from a single c loud. Because of the scaling properties of H 2 pumping 



(jElitzur et all Il989l ). maser clouds with n ~ 10 7 cm 3 and A^h ~ 10 23 cm 



produce near-optimal inversion and detected radiation in cloud-cloud amplifi- 
cation. Proper motion measurements and comparisons of the disk and outflow 
masers offer a most promising means to probe the structure and motion of 
TOR clouds. 



To produce a TOR with the required angular width, clouds must maintain 
a column density above ~ 10 21 cm -2 during their rise to height H ~ r. If a 
cloud expands at the speed of sound c s while rising with velocity v, then its 
size will increase by a factor x ~ 1 + (r/R c )(c s /v). We expect c s /v ~ 10 -2 
since c s is lkms -1 at a temperature of 100 K and the velocity scale of the 
cloud outflow is expected to be comparable to the Keplerian velocity, t>K = 
208 (M. 7 /r pc ) 1 / 2 kms -1 . With R c /r ~ 10 -3 , the 3-D expansion of such a cloud 
during the rise to H ~ r will cause its size to increase by a factor x ~ 10 and 
its column density (oc M c \/R c 2 ) to decrease by a factor x 2 ~ 100, as required. 
Resistance to tidal shearing dictates that the density inside the cloud increase 
as 1/r 3 as it gets closer to the central SBH (eq. [TT]). A cloud of a given density 
can exist only beyond a certain distance from the SBH; only denser clouds can 
survive at smaller radii, leading to a radially stratified structure. Similarly, the 
Keplerian velocity increases as 1/r 1 / 2 as the SBH is approached. Typical BLR 
cloud densities and velocities, discussed elsewhere in these proceedings by H. 
Netzer, occur at r ~ 10 16 -10 17 cm. It is possible that the BLR inner boundary 
occurs where the clouds can no longer overcome the SBH tidal shearing. 
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5.2 TOR Outflow — Properties and Consequences 



The kinetic luminosity of the TOR outflow is L™ R = / ^M c \v%i • uqv & dA. The 
integration is over the disk area where TOR clouds are injected out of the plane 
with vertical velocity component v c \. This injection velocity can be estimated 
from the dispersion velocity of molecular material in the disk — an outflow 
can be considered to have been established when the velocity of the ordered 
motion becomes comparable to that of the local rando m motions. From the 



maser observations in NGC 3079 lKondratko et al.l (120051 ) find that the velocity 



dispersion in a small region (< 5xl0 16 cm) of strong emission is ~14 kms^ 1 , 
which can be considered typical of the local random motions. Therefore, we 
parameterize the injection velocity as v c \(r) = 10 kms" 1 xvQ(Rd)u(r), where 
veiRd) = i>ci(-Rd)/10kms _1 is expected to be of order unity and where u = 
v(r)/v(Rd). Without a full solution of the clumpy wind problem, which is far 
from becoming available, the dimensionless profile u(r) of velocity variation 
across the disk plane remains unknown. However, it is quite certain that u(r) 
decreases with r from its maximum of u(Rd) = 1; in particular, if the turbulent 
velocity is a fraction of the local Keplerian velocity then u = (Rd/r) 1 ^ 2 . With 
M cl nc= m p Ni{N c and Rd from eq. [2] and taking T sub = 1500 K, the outflow 
kinetic luminosity is 

L™R = 7xl0 35 LH 2 N^ 23 vl(Rd)xI 3 ergs" 1 . (12) 



Here N^\ 3 = MqNh^ is the column density through all clouds in the outflow 
equatorial plane in units of 10 23 cm -2 and I3 = J Y (u 3 /y)dyxY/(Y — 1) is an 
unknown dimensionless factor of order < 1. This result shows that the kinetic 
luminosity of the TOR cloud outflow is generally negligible in the overall 
energy budget of the AGN. This estimate for L^ OR does not lend support to 
the suggestion known as the "AGN-Galaxy connection" that outflows from the 
AGN may significantly affect the host galaxy. However, although the energy 
outflow is not a significant factor, the mass outflow is. A similar calculation 
for the mass outflow rate of TOR clouds yields 

M G T U ° R = 0.02 L% 2 A^ 23 v 6 (R d ) x h M yr" 1 (13) 



where I\ = f x (u/y)dyxY/(Y — 1) is another unknown factor of order < 1. 
Since the clouds are likely imbedded in a continuous wind, the actual TOR 
outflow rate is probably somewhat higher. Steady-state mass conservation 
requires accretion into the TOR with a rate M™ R > M™ R ; if this con- 
straint is not met, the TO R mass is depleted within a few Keplerian orbits 
(jElitzur fc Shlosmanl . 120061 1 . Mass accreted into the TOR and not injected 
into the outflow finds its way into the black hole, and is the fuel for the AGN 
luminosity. Therefore M™ R = + M™, where is the black hole 
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accretion rate, so = aM™ where a < 1. In addition, M™ is related to 
the AGN bolometric luminosity via L = rjM^c 2 , where 77 ~ 0.1 is the grav- 
itational conversion efficiency. We end up with L = eM™ R c 2 , where e = rja, 
and the accretion into the TOR and the outflow away from it are related via 



M TOR e 

TTfk = 717^3-6 (R d )xh- (14) 



This ratio cannot exceed unity, yet it increases when L decreases with all 
other factors remaining constant. Therefore, when the luminosity decreases 
the bound M™ R < Mj° R is eventually violated, implying that the system 
cannot sustain the cloud outflow rate inferred from the TOR properties. The 
wind injection velocity v c i(R d ) is likely to change too when L decreases be- 
cause the dust sublimation radius moves closer to the black hole (eq. EJ), where 
all velocities become larger. As with any location in the accretion disk, the 
velocity scale of dynamic motions is determined by the local Keplerian ve- 
locity v K (R d ) = (GM,/R d fl 2 oc (LEdd/L 1 / 2 ) 1 / 2 . If we assume that v cX (R A ) is 
proportional to vk(R&), eq. [TH becomes 



^u? R e ' ( ^Edd\ 1/2 



<23 (15) 



where the modified efficiency factor e' absorbed all additional proportionality 
constants, including Ii. At low accretion rates, the accretio n may sw i tch to 



a rad i atively inefficient mode (RIAF; for recent reviews see iNarayanl . 12002 



M. namely, the efficiency i decreases too and we Veattl^ect the 



ratio e' / L l / A not to vary too much. The TOR equatorial column density must 
remain iVjy 33 > 1 as a definition of the AGN torus, to provide the optical 
depth required by the obscuration and observed IR emission. All in all, de- 
creasing the luminosity at a fixed black-hole mass {L Edd ) causes an increase 
of M™ R /M™ R , and when this ratio exceeds unity the TOR disappears. The 
torus should disappear in low-luminosity AGN. If we take e ~ 0.01 and all 
other parameters ~ 1 in eq. [HI the torus disappearance should occur at bolo- 
metric luminosities below ~ 10 42 ergs -1 . Since the Eddington luminosity of a 
1O 7 M black- hole is 10 45 ergs -1 , the torus disappearance can be expected at 
a typical Eddington ratio L/L Edd ~ 10 -3 . 



5.3 The AGN Low-Luminosity End 



A key prediction of the wind scenario is that the torus disappears at low 
bolometric luminosities (< 10 42 ergs -1 ) because mass accretion can no longer 
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sustain the required cloud outflow rate, i.e., the large column densities. Ob 
servations seem to corroborate this prediction. In an HST stu dy of a complete 



sampl e of low- luminosity (< 10 42 ergs x ) FR I radio galaxies, IChiaberge et al. 



( 119991 ) detected the compact core in 85% of sources. Since the radio selec- 
tion is unbiased with respect to the AGN orientation, FR I sources should 
contain similar numbers of type 1 and type 2 objects, and Chiaberge et 
al suggested that the high detection rate of the central compact core im- 
plies the absence of an obscu ring torus. This suggestion was corroborated by 



Whysong &: Antonuccil (120041 ) who demonstrated the absence of a dusty torus 



in M87, one of the sources in the FR I s ample, by placing stri ngent limits on its 



thermal IR emission. Observations by iPerlman et al.l (120071 ) further solidified 



this conclusion. The IR radiation they detect in the M87 core comes primarily 
from non-thermal jet emission, with only a trace of thermal emission that is 
much weaker than what would be expected from an AGN torus and that can 
be attributed to neighboring dust. The behavior displayed by M87 appears to 
be common in FR I sources. Van der Wolk et al report in these proceedings the 
results of high resolution 12 /j,m imaging observations of the nuclei of 27 radio 
galaxies, performed with the VISIR instrument on the VLT. The reported 
observations provide strong confirmation of the torus disappearance in FR I 
sources. They show that all the FR I objects in the sample lack dusty torus 
thermal emission, although they have non-thermal nuclei. Thermal dust emis- 
sion was detected in about half the FR II nuclei, which generally have higher 
luminosities. In contrast, in almost all broad line radio gala xies in the s ample 



the observations detected the thermal nucleus. Significantly, I Ogle et al.l (120071 ) 
find that most FR I and half of FR II sources have L/L-^d^ < 3 • 10~ 3 , while 
all sources with broad Balmer lines have L/L-^dd > 3 ■ 10~ 3 . 

The low-luminosities source s known as L INER s provide additional evidence 



for the torus disappearance. iMaoz et al.l (120051 ) conducted UV monitoring of 



LINERs with L < 10 42 ergs -1 and detected variability in most of them. This 
demonstrates that the AGN makes a significant contribution to the UV radia- 
tion in each of the monitored sources and that it is relatively unobscured in all 
the observed LINERs, which included both type 1 and type 2. Furthermore, 
the histograms of UV colors of the type 1 and 2 LINERs show an overlap be- 
tween the two populations, with the difference between the histogram peaks 
corresponding to dust obscuration in the type 2 objects of only ~ 1 magnitude 
in the R band. Such toroidal obscuration is minute in comparison with the 
torus obscuration in higher luminosity AGN. The predicted torus disappear- 
ance at low L does not imply that the cloud component of the disk wind is 
abruptly extinguished, only that its outflow rate is less than required by the 
IR emission observed in quasars and high-luminosity Seyferts. When M™ R 
drops below these "standard" torus values, the outflow still provides toroidal 
obscuration as long as its column exceeds ~ 10 21 cm -2 . Indeed, Maoz et al find 
that some LINERs do have obscuration, but much smaller than "standard". 
Line transmission through a low-obscuration torus might also explain the low 
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polarizations of broad Ha lines observed by iBarth et al.l (119991 ) in some low 
luminosity systems. 



If the toroidal obscuration were the only component removed from the system, 
all low luminosity AGN would become type 1 sources. In fact, among the LIN- 
ERs monitored and found to be variable by Maoz et al there were both sources 
with broad Ha wings (type 1) and those without (type 2). Since all objects are 
relatively unobscured, the broad line component is t r uly m issing; in the type 2 
sources in this sample. Similarly, Panessa &: Bassa"n3 ( 2002 ) note that the BLR 
is weak or absent in low luminosity AGN, and Laorl ( 20031 ) presents arguments 
that some "true" type 2 sources, i.e., having no obscured BLR, do exist among 
AGNs with L < 10 42 ergs" 1 . The absence of broad lines in these sources cannot 
be attributed to toroidal obscuration because their X-ray emission is largely 
unobscured. These findings have a simple explanation if when L decreases fur- 
ther beyond the disappearance of the TOR outflow, the suppression of mass 
outflow spreads radially inward from the disk's dusty, molecular region into 
its atomic, ionized zone. Then the torus disappearance, i.e., removal of the 
toroidal obscuration by the dusty wind, would be followed by a diminished 
outflow from the inner io nized zone and di sappearance of the BLR at some 
lower luminosity. Indeed, lOgle et al.l (120071 ) find that all sources with broad 
Balmer lines have L/L-^dd > 3 • 10~ 3 , as noted above. Such an inward pro- 
gression of the outflow turnoff as the accretion rate decreases can be expected 
naturally in the context of disk winds because mass outflow increases with 
the disk area. A diminished supply of accreted mass may suffice to support 
an outflow from the inner parts of the disk but not from the larger area of its 
outer regions. And with further decrease in inflow rate, even the smaller inner 
area cannot sustain the disk outflow. Since the accreted mass cannot be chan- 
neled in full into the central black hole, the system must find another channel 
for release of the excess accreted mass, and the only one remaining is the radio 
jets. Indeed, iHol (120021 ) finds that the AGN radio loudness 1Z = L rSkdio /L opt 
is inversely corr elated with the mass acc r etion rate L/L-Edd, a correlation fur- 
ther verified by I Greene. Ho. fc Ulvestadl (120061 ) . This finding is supported by 



Sikora et al.l (120071 ). who have greatly expanded this correlation and found an 



intriguing result: 1Z indeed increases inversely with L/L Edd , but only so long 
as L/L^d remains > 10 -3 . At smaller accretion rates, which include all FR I 
radio galaxies, the radio loudness saturates and remains constant at TZ ~10 4 . 
This is precisely the expected behavior if as the outflow diminishes, the jets 
are fed an increasingly larger fraction of the accreted mass and finally, once 
the outflow is extinguished, all the inflowing material not funneled into the 
black hole is channeled into the jets, whose feeding thus saturates at a high 
conversion efficiency of accreted mass. It is important to note that radio loud- 
ness reflects the relative contribution of radio to the overall radiative emission; 
a source can be radio loud even at a low level of radio emission if its overall 
luminosity is small, and vice versa. 
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Fig. 12. Conjectured scheme for AGN evolution with decreasing accretion rate. 



The evolutionary scheme just outlines is sketched in figure [121 A similar anti- 
correlation between radio loudness and accretion rate exists also in X-ray 
binaries. These sources display switches between radio quiet states of high 



X-ray emission and radio loud states with low X-ray emission (IFender et al. 



20041 ). While in X-ray binaries this behavior can be followed with time in a 
given source, in AGN it is only manifested statistically. Comparative studies 
of AGN and X-ray binaries seem to be a most useful avenue to pursue. 



5.4 Limitations 



The scheme just outlined was not derived from a detailed dynamical model for 
clumpy disk outflow in AGN. Such models do not yet exist. Instead, the results 
follow from considerations of basic conservation laws within the framework of 
the wind paradigm for the AGN torus. The conclusion that the TOR should 
disappear at low accretion rates is almost unavoidable because the outflow 
rate decreases more slowly than does the luminosity. In eq. [HI our ignorance 
of all dynamics details is lumped into the efficiency factor e. Although this 
factor is largely unknown, the torus disappearance conclusion would be off 
only if e varied with luminosity in a way that completely offset the explicit 
L- dependencies in eqs. HH and LEI Observations indicate that this is not the 



case. The most recent review by iHol (120081 ) presents extensive observational 



evidence for the disappearance of the torus and the BLR in low luminosity 
AGN. 

The lack of explicit dynamics model is both the weakness and strength of 
the presented arguments. Our results provide constraints that have to be met 
by any detailed model constructed to explain the AGN torus within the disk 
wind paradigm. It should be noted that we were able to estimate only the 
properties of the TOR outflow and deduce that this outflow must be turned 



33 



off at fow fuminosities. The BLR disappearance is then deduced from the 
observational evidence that the low luminosity AGN do not be come exc l usivel y 
type 1 objects. A possible dynamical explanation is offered by iNicastrol (120001 ). 
whose work on the role of radiation pressure provides the earliest prediction 
for BLR disappearance at low L/L^d- 



5.5 Final Speculations 



While outflow considerations set limits on the AGN low luminosity end, they 
do not provide any constraints on the TOR at high luminosities. From obser- 
vations, the fraction f 2 of obscured AGN decreases with bolometric luminosity 
(see sec. 13.21) . As is evident from figure [TT| this can be produced by either a 
decrease of a at constant Ao or a decrease of Ao at constant a . The torus 
radial dimensions are set from its inner radius Rd, which increases with L due 
to dust sublimation (eq. [2}. Therefore, the decreasing-a option would arise if 
the TOR height were independent of luminosity or increased more slowly than 
L 1 / 2 , the decreasing-Ao option would arise if the TOR outer radius were inde- 
pendent of luminosity or increased more slowly than L 1 / 2 . The outflow outer 
radius is plausibly confined to within the SBH gravitational influence radius, 
R Q < 35pcx (M.7/^ 2 ) 1 / 3 (sec. 5). As noted in sec. I3.3[ an upper bound on the 

1 /2 

torus outer radius set conservatively from IR observations is R Q < 12 L 4 g pc. 
So whenever L 45 > 9(M. 7 /f2 2 ) 2 / 3 , the TOR outflow is cut at an outer radius 
smaller than that containing the typical number of clouds according to IR 
observations, providing a possible explanation for the decreasing-Ao option. 

Finally, the AGN bolometric luminosity implies accretion into the black hole 
at the rate Af™ ~ 0.1 L 45 M Q yr 1 , assuming radiation efficiency rj ~0.1 (see 
sec. 15.21) . The AGN phase lasts probably ~10 7 -10 8 years, thus the overall mass 
accreted is ~ 10 6 -10 7 L 45 M . Since Seyfert luminosities do not exceed ~10 45 
ergs -1 at most, the entire mass they typically accrete from the host galaxy 
during an AGN episode is no more than that contained in a single giant 
molecular cloud. In contrast, QSO luminosities reach upward of 10 48 ergs" 1 . 
If QSO lifetimes are similar, the mass fed into the SBH is comparable to that 
of an entire galaxy. This suggests that the quasar and Seyfert nuclear activities 
could be triggered differently. Starting a quasar might require interaction with 
another galaxy to trigger accretion of the prodigious amounts of fuel that can 
only be supplied by a full galactic mass. In contrast, a Seyfert episode may 
represent the occasional capture of a giant molecular cloud by the SBH that 
is now known to exist at the center of most galaxies. 
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6 Concluding Remarks 



The clumpy nature of the obscuring torus, long ago deduced from theoretical 
considerations and now being verified by observations, dictates a fundamental 
change in the approach to some basic issues regarding AGN. Foremost among 
them is analysis of the type 1 and type 2 source statistics. Because the viewing 
angle determines the obscuration only in the sense of probability, and not ab- 
solutely, translating source statistics into torus parameters becomes a difficult 
task. In addition to the torus angular width a, which has been used as the 
sole parameter of such analysis, it involves also the mean number of clouds 
A/o, making it impossible to determine both of them from the single mea- 
sured quantity f'2- In particular, the observed decrease of f% with increasing 
luminosity can be attributed to a decrease of either a or A/q or both. Mean- 
ingful constraints on the torus physical properties can be placed only through 
combination of different types of data. Comparisons of UV/optical and X- 
ray properties in individual sources and in samples offer especially promising 
analysis modes for deducing detailed torus properties. 

The low-luminosity end of AGN offers a rich variety of interesting phenomena. 
The consecutive disappearances of the torus and the BLR at low luminosities, 
or low Eddington ratios, and the simultaneous increase in radio loudness fit 
naturally into the outflow paradigm for the AGN torus. Detailed observational 
and theoretical studies of these phenomena would enhance significantly our 
understanding of the black-hole environment. 
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